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The decoherence of two qubits, coupled with one-mode cavity separately, has been investigated 
exactly. The results show that, for the resonant case, the decoherence behavior of system is similar 
to Markovian case when the coupling strength is weak, while the concurrence vanishes in finite time 
and might recover fractional initial entanglement before it permanently vanishes when the coupling 
strength is strong. And for detuning case, the entanglement could periodically recover after a period 

qq ■ of time from its disappearance. These results are quite different from that of system subjected to 

f"^ ' Jaynes-Cummings model. 
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I. INTRODUCTION 

Recently, many works have been devoted to investigating the decoherence behavior of entan gled qubits, invoked by 



coherence behavior ot entangled qubits, invoked by 
the phenomenon, termed as "entanglement sudden death" (ESD) 0, i, i, i, i, i, 0, !, ©, lHI, El 0, El El • It 
is shown that the interaction, between qubits with environment, could lead to finite-time disentanglement, which is 
O ■ quite different from the case of continuous variable two-atom model discussed by Dodd and Halliwelljll. 1 161]. 

The decoherence of two atoms, coupled with one- mode cavity separately, has been investigated in Ref.|3J] with 
Jaynes-Cummings model, which neglects counter rotating terms corresponding the emission and absorption of virtual 
photon without energy conservation, and is widely used in quantum optics (l7l [Lsl. Il9j. Generally, the rotating-wave 
approximation(RWA), which neglecting counter rotating, is justified for small detunings and small ratio of the atom- 
field coupling divided by the atomic transition frequency. In atom-field cavity systems, this ratio is typically of the 
order 10~ 7 ~ 10~ 6 . Recently, cavity systems with very strong couplings have been discussed [20l|. In solid state 
s ! ' systems, the ratio may become so large as to consider the effect of counter-rotating wave terms [2 lj. In this paper, we 
investigate the influence of counter-rotating wave terms, which are neglected in Ref.jj], on the decoherence behavior 
of two atoms coupled with one-mode cavity separately. 

In section[TlJ the reduced non-perturbative quantum master equation of atom has been derived and its exact solution 
is obtained. In section ITTTl the decoherence of two initially entangled qubits has been discussed. The conclusion will 
be given in section HVl 

II. MODEL AND EXACT SOLUTION 
A. non-perturbative master equation and its Exact solution 
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Now we restrict our attention to two noninteracting two- level atoms A and B, which are in a perfect one- mode 
cavity, respectively. First, we consider the subsystem Hamiltonian of one atom coupled to one cavity field mode as 

H = H a + H f + H af (1) 

where 

Ha = W y (2) 

Hf = ujoJ a (3) 
H af = g{a + + a-){at + a) (4) 

where ojq is the atomic transition frequency between the ground state |0) and excited state |1). a z = |1)(1| — |0)(0|, 
c + = |1)(0| and er_ = |0)(1| are pseudo-spin operators of atom, and a are creation and annihilation operators of 
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the cavity field mode corresponding frequency u>. And g is the coupling constant between the transition |1) — |0) and 
the field mode. 

The reduced non-perturbative quantum master equation of atom could be obtained by path integrals [22j 
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Pa = -i£ aPa - / ds{C af e- iCo ^C af e- ij: ^ s -^} f p a 
Jo 



(5) 



where Co, C a and C a f are Liouvillian operators defined as 



and 



/(*) 



t'A 

exp(iSt) — 1 
iS 



(6) 



C p = [H a + H f ,p] 

C a p EE [H a ,p] 
C afP EE [H a f,p] 

and (...) / stands for partial trace of cavity mode. 

If the cavity field is initially in vacuum state , the non-perturbative reduced master equation of the atom could be 
derived from Eq.© 

^ Pa = ~g 2 (a R + f(t)) p a - 2i (loo - 9 2 a' + g 2 f{t)) J aPa 

+g 2 (a + f*(t)) J +P a + g 2 (a* + /(f)) J_p a + 2g 2 a R K +Pa 
+2g 2 (a R - f R (t)) Kop a + 2g 2 f R (t)K_p a 

Where Jo, J+, J-, Kq, K + and K- are superoperators defined as 

7 - \° z 

J()Pa = \~£iPa 

J+Pa = V+Pa&+ 
J- Pa = O-PaO- 

K Pa = (a + a-p a + Pa a + a- - pa)/2 

K+Pa = O+PaO- 
K-p a EE 0_p a <7 + 

1 — exp(-iAt) 



(7) 

(8) 



where A = to + cjo, 6 = u>o — u>. a , a 1 , a* and / , f 1 , f* are real part, image part and conjugate of a and of f(t), 
respectively. 

Using algebraic approach, the formal solution of Eq.© is obtained dHHH 



p a (t) = exp(-T k )fexp 



dt(eoJo + s+J+ + E-J-) 



xTexp 



dt(u K + v+K+ + v-K-)p a (0) 



(9) 



where t is time ordering operator, e — ~2z (ui — g 2 ® 1 + g 2 / 1 ) > £+ 
2g 2 (a R - f R ), v+ = 2g 2 a R , V- = 2g 2 f R , T k = g 2 [a R + F R ) and 



g 2 (a + f*), £_ = g 2 (a* + f), u 



1 — exp(— iAt) — iAt _p r 
adt = —7: ee a + ia 



F = I f(t)dt 
h 



A 2 

1 + iSt — exp(iSt) _ R 



S 2 



F R + iF 1 



(10) 
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where a R , a 1 , a* and F R , F 1 , F* are real part, image part and conjugate of a and of F, respectively. 
Using the decomposition of SU(2) operator, the exact solution of master equation Eq.© is obtained 



P(t) = 



Pa{t) 



m 



_(lpl 1 (Q)+m f %>(0) «Pi°(0) +VP° a 1 (0) 



(11) 
(12) 
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e ~h/2j_ 



gJo/2 + e -io/2j + j_ 



y 



e - 30 /2 J+ 



(13) 
(14) 
(15) 
(16) 



where j+, jo, J- and k+, ko, fc_ satisfy the following equation [l£ 

X+ = /!+ - ll-X\ 

X = no - 2/i_X + 



p X A 



(17) 
(18) 
(19) 



H = £ for X = j and p, = v for X = fe. 



Concurrence 



Throughout the paper, we use Wootters concurrence |24| . For simplicity, we assume that the two subsystems have 
the same parameters. The concurrence of the whole system could be obtained[6J] 



Q = max{0,ci,c 2 },(f = 
ci = 2e~ 2rfc (^/p 2 3P32 - VP11P44) 



(20) 



c 2 = 2e 



-2r t 



(V Pupa - VP22P33) 



corresponding to the initial states of |$) = /3|01) + 77! 10) and |^) = /3|00) + 771 1 1} , respectively. Where j3 is real and 
< /3 < 1, 77 = 1 77 1 and f3 2 + \r/\ 2 = 1. The reduced joint density matrix of the two atoms, in the standard product 
basis B = {|1) = |11), |2) = |10), |3) = |01), |4) = |00)}, could be written asQ 



-2r fc 



/Pll P14 

/0 2 2 P23 

P32 P33 

V P41 p 4 4 



(21) 



here the diagonal elements are 



P11 
P22 

P33 
P44 



l 2 pn{0) + lmp 22 {0)~ 
lppu(0) + lmp 22 (0) - 
lppn(0) + pmp 22 (0) 
p 2 pn(0) +pnp 22 (0) - 



mlp 33 (0) + m 2 p 44 (0) 
mpp 33 {0) + mnp Ai (0) 
- nlp 33 (0) + nmpufi) 
■ npp 33 (0) + n 2 p 44 (0) 



and the nondiagonal elements are 



pu 
P23 

P32 
PA1 



x 2 p 14 (0) + xyp 23 (0) + yxp 32 (Q) + y 2 p 4 i(0) 
xrp 14 (0) + xqp 23 (0) + yrp 32 (0) + yqp^O) 
rxpu(0) + ryp 23 (0) + qxp 32 (0) + qyp 4 i{0) 
r 2 pi4 (0) + rqp 23 (0) + qrp 32 (0) + q 2 p 41 (0) 
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III. NUMERICAL RESULTS AND DISCUSSION 



In order to compare the results with that of two Jaynes-Cummings atoms in Ref. |3| , we primarily investigate the 
decoherence for resonant case u> — u>q. 

First, we focus on the decoherence of two qubits with initial state of |<I>). For the RWA model in Ref.Q, FigQ]shows 
that the concurrence periodically vanishes and revives. And the change of C$ against 1 is symmetrical because of 
the symmetry of the initial state |$). The decoherence, of the non-RWA model in this paper, has been shown in Fig|2] 

(A) As wo = 1.5(7, FiglUa) reveals that the concurrence C$ decreases to zero in a finite time, vanishes for a period of 
time, revives with small amplitude and then vanishes permanently. This characteristic will hold on for more stronger 
coupling constant. 

(B) From FigJ^b), we could find that, when lu = 3g, the concurrence Cj> first decreases to a certain value and 
maintain it for a period of time before it vanishes, like a rumple, then it revives with a smaller amplitude than that 
in Fig[2ja) and vanishes permanently. As the coupling constant g decreases, there are more rumples and a smaller 
amplitude revival. 

(C) When ujq = 30g, Figl2tc) exhibits that the concurrence C$ decreases monotonically and exponentially to zero 
without revival of entanglement. This characteristic will hold on for more weaker coupling constant. 

Then, we focus on the decoherence of two qubits with initial state of For the RWA model in Ref. [l2| . Figj3]shows 
that the entanglement represented by can fall abruptly to zero, and will remain zero for a period of time before 
it recovers. The length of time interval for the zero entanglement is dependent on the degree of initial entanglement. 
The time interval for /3 2 < 1/2 is longer than that for /3 2 > 1/2. There are periodical disappearance and revival of 
entanglement in time scale. For the non-RWA model, the decoherence of the system was also categorized into three 
cases 

(A) When ujq = 1.5g, FigfJJa) reveals that, for 1 > 1/2, the concurrence decreases exponentially to zero, 
remain zero for a period of time, revives fractional initial entanglement and then vanishes permanently, while 
vanishes permanently after a finite time for (3 2 < 1/2, similar to the Markovian case[l[. 

(B) As uiq = 3g, Figdjb) shows that the entanglement represented by dp has a similar behavior to that of C$ for 
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1 > 1/2 in FigJ^b). In contrast, for /3 2 < 1/2, it decreases smoothly to zero at finite time, then vanishes permanently 
after a small amplitude entanglement revival link that for 1 > 1/2. 

(C) When ojq = 30<?, Fig|4jc) exhibits that it first decreases to zero at short time, then vanishes permanently 
after a small amplitude entanglement revival. Unlike the two cases above, the evolution behavior of concurrence C-^ 
becomes symmetric, like that of C$, because of the strong interaction of atom with reservoir through the emission 
and absorption of virtual photon. 

From Fig|2] and Fig[4] we find that the entanglement will decreases to zero finally for no-RWA and there are no 
periodical disappearance and revival of entanglement, like that in FigfTJand FigJ3]for RWA case. It is the enhancement 
of spontaneous emission, as an atom resonantly coupled with a cavity, leads to the disappearance of entanglement (25[. 
We could conclude that the neglect of counter-rotating wave terms in Hamiltonian of RWA model leads to the different 
characteristics of C<s> . 

Finally, we investigate whether there are periodical disappearance and revival of entanglement for detuning case. 
Concurrence C$ and for detuning case as a function of gt and (3 2 are shown in Fig[5] and Fig[6j respectively. 

From FigJS] and FigJSl we could find that the initial entanglement will recover after a period of time from its 
disappearance because the spontaneous emission could be greatly inhibited, as an atom non-resonantly coupled with 
cavity mode[25|. The revival time interval is dependent on the detuning S. The bigger the detuning is, the shorter 
the time interval is. We also find that the entanglement will change little when the detuning is large and coupling 
strength is weak. 



IV. CONCLUSION 



The decoherence of two initially entangled atoms, coupled with two one-mode cavities separately, has been discussed 
exactly. The results show that the decoherence behavior of two qubits is dependent on the coupling strength and the 





01 2345678 

gt 

FIG. 6: Concurrence C* and its contour for non-resonant case as a function of gt and 1 with ujo = 10<;,<5 = O.ltJo- 

detuning between atom transition frequency and the cavity mode. 

Firstly, there are no periodical disappearance and revival of entanglement for resonant case like that for Jaynes- 
Cummings model in Ref. [|[. Secondly, for detuning case, the entanglement could periodically recover after a period of 
time from its disappearance. Thirdly, for resonant case, the decoherence behavior of system is similar to Markovian 
case when the coupling strength is weak, while the concurrence will vanishes in finite time and might recover fractional 
initial entanglement before it permanently vanishes when the coupling strength is strong. 

The results also exhibit that the RWA in Hamiltonian leads to the existence of revival of entanglement for resonant 
case and it might be improper to take RWA when the interaction between atom and external field is correlated and 
strong. 
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